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Long distance coupling of resonant exchange qubits 


Maximilian Russ and Guido Burkard 
Department of Physics, University of Konstanz, D-78457 Konstanz, Germany 

We investigate the effectiveness of a microwave cavity as a mediator of interactions between two 
resonant exchange (RX) qubits in semiconductor quantum dots (QDs) over long distances, limited 
only by the extension of the cavity. Our interaction model includes the orthonormalized Wannier 
orbitals constructed from Fock-Darwin states under the assumption of a harmonic QD confinement 
potential. We calculate the qubit-cavity coupling strength in a Jaynes Cummings Hamiltonian, and 
find that dipole transitions between two states with an asymmetric charge configuration constitute 
the relevant RX qubit-cavity coupling mechanism. The effective coupling between two RX qubits in 
a shared cavity yields a universal two-qubit iSWAP-gate with gate times on the order of nanoseconds 
over distances on the order of up to a millimeter. 


I. INTRODUCTION 

Quantum computation with single-electron spins con¬ 
fined in semiconductor quantum dots (QDs^has been in¬ 
vestigated within a wide range of implementations yield¬ 
ing long decoherence times on the order of ^ /rs.^HS Gal¬ 
lium Arsenide (GaAsj^ and silicon (Si^ are the most 
common choices of host materials for the QDs. By com¬ 
parison, charge qubits decohere much faster (within ^ 
npni) 

due to the strong coupling between the charge and 
the electromagnetic fields in the environment.^ There¬ 
fore, the aim of many implementations is a high protec¬ 
tion against electrical noise to achieve long qubit decoher¬ 
ence times in order to be able to perform as many qubit 
operations (gates) as possible while the qubit is coherent. 
Experiments show high-fidelity state manipulation and 
long coherence times for scalable implementations using 
single or multiple QDs.li^l^^ One promising candidate is 
the resonant exchange (RX) spin qubit, a modification 
of the exchange-only qubit which allows for all electric 
control of the qubit for the price of a triple quantum dot 
(TQD) scheme. This implementation possesses a high 
robustness against electrical noise at the sweet spot ,221211 
but is still susceptible to electromagnetic fields at the 
resonance frequency allowing for additional q ubit c ontrol 
through radio-frequency or microwave signals.fi 222 l 
shown that two-qubit gates between RX qubits can be 
implemented with the exchange coupling, using a single 
exchange pulse. 

The coherent transport of quantum information, e.g., 
long distance entanglemenpSi, combined with state 
preparation and read-out is investigated within a wide 
range of implementation^^ in quantum optics. Adapting 
techniques from this field, a long distance coupling be¬ 
tween two solid-state QD spin-qubits can be envisioned, 
which could complement the existing short-range interac¬ 
tions, such as the exchange coupling between nearby QD 
spin-qubits. In order to achieve long distance coupling 
between spin-qubits, a long-range interaction is needed, 
e.g., the coupling of the qubits to an electromagnetic field 
with specific photon mode^^ or by the coupling of the 
qubit to a ferromagnet.^In this paper, we focus on long¬ 
distance coupling between RX qubits based on the elec¬ 



FIG. 1. Schematic illustration of the proposed setup for long 
distance interaction. The setup consists of two linearly ar¬ 
ranged TQDs (RX qubits) inside a superconducting strip-line 
cavity. For an optimal setup the two qubits should be lo¬ 
cated at the held maxima (anti-nodes) in order to achieve 
a strong qubit-cavity coupling. The green environment sur¬ 
rounding the RX qubits illustrates the coupling of the qubit 
to the electromagnetic held of the cavity (blue). As a result, 
the RX qubit can be coupled by using photons as mediators. 
The gray layer illustrates the two-dimensional electron gas 
(2DEG) in which the quantum dots are embedded. 


tromagnetic fields in a microwave cavity (Fig. [^. 

Goupling spin-qubits via electromagnetic cavities has 
been proposed for spin-qubits in nitrogen vacancy (NV) 
centers in diamond. NV centers can be coupled with 
photons in the optical spectrum,l^2Ell while QD spin- 
qubits usually react to microwave or radio-frequency 
signals. Since microwave cavities with a high finesse 
exist ,1^ this long-range coupling is feasible between QD 
spin-qubits in solid-state materials, e.g., coupling spin 
qubits to photons in a cavity via electric dipole or gate 
potentials.23311351 Experiments have shown evidence of a 
strong coupling between qubits in single or double QDs 
to a micr owave cavity through the charg^SlHlIl qj- 
spin 140142 Here, we present an implementation of such 
a long distance interaction between two RX qubits en¬ 
abled in a TQD scheme at time scales on the order of 
nanoseconds. 

This paper is organized as follows. In Section [ll] we 
introduce our model for the long distance interaction. 
Subsequently, in Section [m] we calculate the associated 
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transition dipole matrix elements used to determine the 
corresponding qubit-cavity coupling parameter and iden¬ 
tify the underlying coupling mechanism. Finally in Sec¬ 
tion m we combine these results to present a step-by- 
step prescription for a universal two-qubit iSWAP-gate be¬ 
tween two RX qubits in a shared cavity. We conclude in 
Section [V| with a summary and an outlook. 

II. MODEL 

We consider two linearly arranged triple QDs (TQDs), 
where each QD has a single available orbital, occupied 
by three electrons. Both TQDs are assumed to lie in a 
superconducting microwave cavity with a single available 
photon mode (see Fig. [^. We describe this system with 
the Hamiltonian 

2 

H = ^(iJ,+7Ji„t.*)+Tfcav, (1) 

z=l 



where Hi describes the dynamics of the electrons in the 
i-th isolated TQD, is the interaction between the 

electrons in the i-th TQD and the photons in the cavity 
and Hcav describes the photons in the cavit y. 

For the qubit we use the RX HamiltoniarP2l^ Hi de¬ 
rived from the three-site extended Hubbard Hamiltonian 
which describes a linearly arranged TQD (see Fig.[^(a)). 
We work in the RX regime in which only the charge 
states (1,1,1), (2,0,1) and (1,0,2) are accessible. The spin 
qubit lies in the subspace with spin quantum numbers 
S = Sz = spanned by the states 
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( 2 ci,t''W, 

|5'i,i/2) = |vac), 

|5'3,i/2) = |vac), 


uV.tV.t) ’ 
( 2 ) 


where |vac) denotes the vacuum state and cl ^ (ci,o-) cre¬ 
ates (annihilates) an electron in QD i with spin cr. We 
have further used the notations l^) = (|ti) — |4.t))/2, 
|r+) = Itt), and |To) = (|ti) + |4.t))/2 for the singlet and 
two of the triplet states of two electrons. Here, {m,n,l) 
denote the number of electrons in the first (n), second 
(to) and third {1) Q D. In this basis, we obtain for the 

Hubbard HamiltoniarPsHll 


H = 


/ 0 0 ti/2 tr/2 \ 

0 0 V3ti/2 -V3trl2 

ti/2 Vsti/2 A + e 0 

V tr/2 -V^tr/2 0 A-e J 


( 3 ) 


where is the hopping between the left (right) and the 
center QD, e is the energy difference between the outer 


FIG. 2. (a) Schematic illustration of the triple quantum dot 
(TQD) confinement potential V{x). The hopping matrix ele¬ 
ment between QD 1 (QD 3) and QD 2 is denoted by ti {tr). 
We also show the response of the system to two electrically 
controlled bias parameters, the difference e between the en¬ 
ergy levels of the outer two QDs and, the effective difference 
A of the energy levels in QD 2 and the mean of the outer 
QDs (including Coulomb repulsion). Here, we assume that 
the TQD is filled with three electrons. We also include the 
Coulomb repulsion in the center QD, Ec = U — 2Uc- (b) 
Schematic illustration of the orbital wave functions of the 
electrons in the TQD. The inter-dot distances a; (a^) between 
QD 1 (QD 3) and QD 2 need to be sufficiently small to al¬ 
low for a sizeable overlap Si {Sr) between the orbital wave 
functions. 


QDs and A is the effective energy difference between 
the center QD and the outer QDs in which Coulomb re¬ 
pulsion is included (see Fig. (a)). In the RX regime 
(|e| < |A| and ti^r |A ± ej) the states |0) and |1) 
are nearly eigenstates, while the states |2) with a charge 
configuration (2,0,1) and |3) with a charge configuration 
(1, 0, 2) are only virtually occupied and can be eliminated 
via a second order Schreiffer-Wolff (SW) transformation 
yielding a Heisenberg model i?Heis = JiSi- S 2 + JrS^ ■ Ss 
with the exchange energies Ji = tf/{A + e) and Jr = 
tr/{A — e). In its eigenbasis which we consider as our 
logical qubit space, the Hamiltonian takes the form, 

h 

Hrx = (4) 

in which fiwRx = is the res¬ 

onance frequency of the qubit. Periodic driving of the 
detuning parameters e and A allows Rabi transitions be¬ 
tween the two eigenstates that are, hence, interesting for 
cavity quantum electrodynamics (cQED). 

For the resonator we consider a superconducting strip- 
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line cavity whose resonance frequency is in the GHz- 
regim e and matches the energy splitting of the RX 
qubit [milllll Since the relevant dynamics of the cav¬ 
ity mostly depends on the dynamic of a single electro¬ 
magnetic mode near the qubit resonance frequency, the 
Hamiltonian is that of a single mode,^ 

Hcax — -j- —^ , (5) 

where (a) is the bosonic creation (annihilation) op¬ 
erator of a photon with frequency Wph- The associ¬ 
ated eigenenergies are ifcav = ^ph (jT-ph + 1/2), where 
riph = (a^a) counts the number of photons with the cav¬ 
ity frequency Wph- 

For the qubit-cavity interaction, we consider the min¬ 
imal coupling Hamiltonian in the dipole approximation 
near the resonanc^^ 


Hint 


e / h 
meff V2eoeRwph 


1/2 

ep-p{a + al) , 


( 6 ) 


where e (eg) describes the dielectric constant of the ma¬ 
terial (vacuum), V is the volume of the cavity, and 
is the polarization of the photons. In the logical qubit 
subspace, the dipole matrix element of interest is 


9r 


e / h 
2m \2eoeRwph 


1/2 

(0|ep -pII) 


( 7 ) 


which describes the photon-induced transition between 
our two qubit states |0) and |1). 


III. QUBIT-CAVITY COUPLING 
A. Phenomenological approach 

To describe the coupling of a qubit to the cavity, we 
quantize the two detuning parameters e, A. This yields 
e —>■ e + K (a -I- a^) and A —> A + k' (a -I- a^), where k and 
k' are parameters that include both, the amplitude seen 
by the qubit for a given polarization and the vacuum 
amplitude of the electric field, and (a) is again the 
creation (annihilation) operator of a photon in the cav¬ 
ity with frequency Wph- Assuming ^(a-l-al) <C e and 
k' (^a + a'') <C A we can expand the RX Hamiltonian in 
terms of k (a -I- a^) and k' {a + a^). Asa result, we fincP^ 

Hrx = ^[{f>^RX ++ Suj^a^^], ( 8 ) 

represented in the eigenbasis of the unperturbed sys¬ 
tem with the longitudinal coupling Suj^ = —{JSJ + 
3jAj) (a -I- al) /ujrx and the transversal coupling = 
\/3{JSj — jSJ) [a + a’') /wrx- Here, we used J = 
{Jl + Jr)l2., 6J = deJ K + OaJk', j = {Jl — 

Sj = dgj K + OaJ k' and the derivatives d^J = OaJ = 
{eJ - Aj)/{A‘^ d^j = OaJ = {ej - AJ)/{A^ - e^). 



FIG. 3. Two possible arrangements of the electric field E 
(blue arrow) in a cavity, (a) The TQD is arranged parallel to 
the electric field resulting in a finite coupling k to the param¬ 
eter £ and in a vanishing transversal qubit-cavity coupling k' 
to A, due to the long wave length of the cavity photons, com¬ 
pared to ai^r- (b) The opposite case in which the two outer 
gates of the QDs are connected to the same potential and the 
electric field E is aligned from the center QD which results 
in vanishing qubit-cavity coupling k due to a static e and a 
finite k! . 


Neglecting higher-order terms of the expansion which 
correspond to two-photon processes and higher, we ob¬ 
tain, 

77jc = ^0-2 + 5s <^x (a + a^) + {a^a + 1/2) (9) 
with the coupling parameter 

5s = V3{JSj - jSJ)/ujRx (10) 


and a photon-dependent resonance frequency tUg = wrx~ 
{JSJ+3jSj) (a -I- ol) /ujRx- For certain alignments of the 
TQD and the cavity, the coupling parameters can be set, 
e.g., for an alignment as in Fig. ^ (a) we expect k' « 0 
due to the long wavelengths of the cavity photons while 
for the alignments in Fig. (b) k is negligible. In the 
case k' = 0, we obtain 


5s 


= VSn 


.2.2 


(A2 — £2) 

and in the case k = 0, we find 


5e 


= VSn 


.2.2 


(A^ — e^) 


( 11 ) 


( 12 ) 


Considering incoherent and broad-band electromag¬ 
netic fields, we find the same expression Eq. ([^ for the 
RX qubit under the influence of charge noise This 
is not surprising, since in both cases the RX qubit is 
disturbed by electromagnetic helds. Hence, increasing 
the qubit-cavity coupling also increases the coupling to 
charge noise and moving to a sw eet spot, where the qubit 
is robust against this noise,l221^ we expect a weak qubit- 
cavity coupling as a trade-off. This phenomenological 
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model does not provide a microscopic description of the 
coupling parameters k and k' . We consider a more realis¬ 
tic model for an alignment as in Fig.|3|(a) in the following 
Section p^II B[ which will also allow us to estimate k. 

B. Microscopic theory 

Coupling the qubit states of the RX qubit in the TQD 
to the cavity via the emission and absorption of a cav¬ 
ity photon requires a strong electric-dipole transition el¬ 
ement gr between the two qubit states. We consider an 
alignment of the TQDs as shown in Fig. (a) with the 
electric field in the cavity pointing in cc-direction. To find 
a finite electric-dipole transition in this system, the fol¬ 
lowing conditions have to be fulfilled. Firstly, the TQD 
has to be coupled by inter-dot exchange interactions, i.e., 
hopping between neighboring QDs. We find that the 
matrix element is approximately proportional to the en¬ 
ergy splitting between the qubit states (~ titr) which 
matches with past calculations with a double quantum 
dot (DQD).I^ However, there are still two independent 
symmetries which have to be broken for a non-vanishing 
matrix element. 

The first symmetry arises from the spin-conserving na¬ 
ture of the electric-dipole transitions and corresponds 
with inversion symmetry. To distinguish the product 
states of three electrons in a linearly arranged TQD we 
use three quantum numbers.!^ The first two are the to¬ 
tal spin S and its z-component Sz which are identical for 
the qubit states, while the third quantum number is the 
total spin Sq of the two electrons in the outer QDs 1 and 
3 which distinguishes them, iSo | 0 ) =0 and S'o |1) = Ij 
hence, they cannot be transferred into each other by an 
interaction O with [iS'o,0] = 0. The influence of the 
states 12) and |3) with asymmetric charge configuration 
breaks this symmetry for Ji ^ Jr, e.g., through a gate 
detuning e y^Q. 

The second symmetry that needs to be broken is an or¬ 
bital symmetry between the electrons in the outer QDs. 
The orbital wave functions of an electron in QD 1 and 3 
have identical parity under the assumption that all three 
QDs have the same lowest orbital confinement potential. 
To overcome this symmetry either the confinement po¬ 
tentials of the two outer QDs need to be different or the 
inter-dot distance between the TQDs has to be asymmet¬ 
ric, hence, o/ ^ Qr where ajyr) is the inter-dot distance 
between QD 1 (QD 3) and the center QD (see Fig.[^(b)). 

1. Wave functions 

To calculate the transition dipole matrix element 
Qr oc ( 0 | ep • p\\) between the qubit states, we need 
orbital wave functions which are pairwise orthogonal, 
since the isolated orbital wave functions have a fi¬ 
nite overlap (see Fig. (b)). For orthogonal states, 
these finite overlaps Si = {ipi\g) 2 ), Sr = {‘f?,\^ 2 ) and 


Si 3 = ((/3i|<P 3) have to be zero which is not the case 
here. Therefore, we transform the non-orthogonal basis 
{|(^i), |(^ 2 )) |‘^ 3 )} into an orthonormal basis of the Wan- 
nier orbitals {|4>i), |<i) 2 ), |$ 3 )} which fulfill the orthonor¬ 
mality condition (<l>i| = 6 ij with i,j G {1, 2,3}. Here, 

these states are chosen such that they describe the dy¬ 
namics in the i-th QD and converge for long distances to 
the isolated electron wave functions, lim |$i) = \(pi). 

ai,ar —^oo 

Hence, we obtain in the general case, 

l^i) = ^ (l</5i) + ai |v?2) + h |(/33)) 

1 ^ 2 ) = ^(|</52)+a2|V5l)+&2|‘/?3)) (13) 

1 ^ 3 ) = ^ (I</53) + 03 W2) + h l‘/?i)), 

where and bi are bounded, real parameters. However, 
these Wannier states are not uniquely defined, since the 
orthogonality condition yields only a linear equation sys¬ 
tem with three equations, but, with nine independent 
parameters. Three parameters {Ni, N 2 , N^) can be elim¬ 
inated immediately by the normalization condition, 

Ni = ^\ 2ai Si + 26i 513 -I- 2ai bi Sr 

N 2 = \ ‘2,(12 Si + 262 Sr 2(22 ^2 S 13 -|- ^2 "b & 2 i (^^) 

N 3 = \J^ 2a3 Sr 263 S '13 + 203 63 S'; -|- 03 + 63 . 

In a long distance approximation, we neglect the overlap 
between QD 1 and QD 3 and set 513 = 61 = 63 = 0. 
This reduces the number of free parameters by two and 
yields simpl er expressions fo r the normalization parame¬ 
ter, 1 -t- ‘ 2 ci\ Si -\- , A ^3 = -f 2(^3 Sr O 3 and 

N 2 = 1/1 + 202 Si + 2&2 Sr + a\-\- h\. The last param¬ 
eter we adapte from the condition of maximally local¬ 
ized Wannier orbitals where we minimize the localiza¬ 
tion functional^ ^ 

(see Appendix |B 1[ ). We found at lowest order a sim¬ 
ple behavior for the parameters Oi and 03 , in particular, 
ai = iaSi + 0{Sr,S^) and ag = iaSr + 0{Si,S'^) and 
Ca = which yields the final condition 01/03 = Si/Sr- 
As a result, we obtain analytical expressions for the 
parameters, Oi = —25'/, 02 = 5//(l — 2 5/ — 2 5/), 
&2 = 5r/(l — 2 5/ — 2 5/) and 03 = —2Sr within the 
scope of the approximation. A detailed parameter dis¬ 
cussion can be found in Appendix |B 2[ 

2. Transition dipole matrix elements 

Having found pair-wise orthonormal wave functions in 
position space we first calculate the matrix elements of 
the position operator x in the one-particle basis which 
we use to calculate (0|pa;|l) and finally gr- Therefore, 
we express the position operator x in the orthonormal 
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basis of the Wannier functions , |$ 2 ), I'ha)} each 

combined with one of the two spin states H) , |t) in order 
to describe the electron spin dynamics. As a result, we 
obtain 

3 

^ (15) 

i,j=l cr=tt 

where Xij = with i,j G {1,2,3} and the op¬ 

erator c[ ^ (ciy.) in Eq. © creates (annihilates) an elec¬ 
tron in the orthonormalized Wannier orbital with spin 
a G {til- In the next step, we construct the qubit states 
|0), |1) and the asymmetric states |2), |3) in terms of the 
spin Wannier states = l^hi) I o’). Keeping this in 

mind, we now express the position operator x in the 


{|0). |1) > |2), |3)} basis as 





/ 3 

I] X^^ 

0 



\ 

X = 

0 

3 

2=1 

\Jlx 12 






2xii -b X 33 

-a;3i 



V ^^23 

~ Y^^23 

-X13 

Xii -b 2 x 33 

/ 


(16) 


We obtain the elements of the momentum operator in 
the same basis through the relation p = [i?Hub, x\^ 
Here, i?Hub is the full Hamiltonian in this basis given in 
Eq. (§ and the square brackets denote the commutator. 
An analytical expression can be obtained and is shown in 
Appendix |D 2 [ The simplified expression in Eq. (D3) can 
be obtained considering real matrix elements, x,, = x,,- 
with i,j G {1,2,3}. 


3. Interaction Hamiltonian 


We now have almost all the tools for calculating the 
qubit-cavity coupling in iJint- However, due to the 
Schrieffer-Wolff transformation the qubit states in the 
expression for the RX Hamiltonian in Eq. Q are not ex¬ 
actly the states defined in Eq. ([^, in particular they have 
a small contribution from the asymmetric states | 2 ) and 
13). Therefore we have to shift into the qubit basis by 
the same transformation, hence, compute p = e^pe ~^ » 
p — [p, , where S is the SW transformation matrix.1221241 

Considering only the matrix elements in the logical sub¬ 
space of the RX qubit {|0), |1)} and neglecting all other 
elements, we obtain as a result (see Appendix [Pj) 


Px =- 


•\/3m titrS 


h 

titrS 


m 


-Re (xia) (Ty 
dm (xia) 


h L — £2 
V^ty.Im(x23) - V^t;lm(xi2) CTz 


(17) 


photon induced 
transition 



FIG. 4. Schematic illustration of the qubit-cavity coupling 
mechanism. Photons induce dipole transitions between the 
virtually coupled states |2) and |3) (top line) which have an 
asymmetric charge configuration due to the double occupa¬ 
tion of one QD. These two asymmetric states couple weakly 
through some hopping matrix elements, here denoted as ti^r, 
to the qubit states. The arrows in |0) and |1) show only one 
of the possible spin configurations. 


where irrelevant terms proportional to the identity oper¬ 
ator in the qubit space have been omitted. 

Comparing this result for px including non-vanishing 
matrix elements with the one for the bare momentum 


operator in the top-left corner of Eq. (D2), one finds 


that in the latter, only the imaginary part of the single 
electron matrix elements appears, e.g., Im(x 23 ), whereas 
Eq. ( [T7| ) also contains the real part of X 13 . Hence, we 
conclude that the physical mechanism of the coupling 
between the qubit states of the RX qubit is due to the 
coupling to the asymmetric states |2) and |3), since in our 
case Xij G K, with i,j G {1,2,3} (see Appendix [^. In 
particular, the electromagnetic field of the cavity induces 
electronic transitions, where one electron from a doubly 
occupied QD is transferred to the singly occupied QD 
and not the empty one, hence, between the states | 2 ) 
and |3) (see Fig. 0. In this case, we hnd 


Px =- 


\/3r 


titrS 

A 2 - £2^13'^!/- 


(18) 


It should be noted at this point that a non-perfect align¬ 
ment of the TQD or electric field in x-direction, e.g. 
a tilting angle in the xj/-plane, gives rise to an imag¬ 
inary contribution to the dipole transition matrix ele¬ 
ments, hence, allowing for other electronic transitions 
(see Eq. However, this contribution is small if the 

projection on the j/-axis, i.e., the tilt angle, is small. 

As a last step, we directly compute the single particle 
matrix elements in order to find the qubit-cavity coupling 
Pr- The calculation of the associated matrix element X 13 
with the help of the orthonormalized Wannier functions 
|$i) yields 


\/6m 

n 


[tilm(xi 2 ) -b tr-Im(x23)] CTx, 


Px = — 


'/3 m titr-e 


h A 2 — £2 


^^rel /l ^r) SlSj-(Ty 


(19) 
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FIG. 5. Qubit-cavity coupling Qr as a function of the inter-dot 
distance a; for (top) GaAs and (bottom) Si. The parameters 
are chosen as ar = l.la;, e = 0.2meV, A = O.SmeV, B = 
310 mT, LOO = 3.1 fJ-eV and tUph = curx ~ 19.54 /reV^. We use 
hopping parameters ti and ty. from numerical calculation^^ 
for the black dots, while we consider the expression for a DQD 
ti^r — hujQo [(ai,r/ao)^ + (tiJo/wQD)] 5!,r/2 (l — Sf^r) for the 
gray curve.^ Here, oq = -y/li/mwo is the confinement radius 
with Hloqu being the confinement energy of the QDs. Typical 
values are for a GaAs QD cjQD.GaAa = 3.1 meV and for a 
silicon QD cjqd.si = 5.9 meV due to their different effective 
masses, moaAs = 0.067mo and msi = 0.191mo, with mo being 
the free electron mass. For the volume V = Ld\d 2 in the 
vacuum amplitude of the electric field we set the length L = 
1 mm, the width di = 1 /rm and the height d 2 = 100 nm of the 
stripline cavity. The red, blue and green dashed lines indicate 
three appropriate values for quantum gate operations. 


where Orel is close to one. An explicit calculation can be 
found in Appendix Substituting the resulting momen¬ 
tum operator into Eq. §) we obtain for the interaction 
Hamiltonian 


Hint = ffr cTy (a + a'i) (20) 

with the effective coupling strength 


and the vacuum amplitude of the electromagnetic field 
in the cavity Eq = (hujph/2€oeV) ' . Identifying gr = gs 
from Eq. (11), we find 


ai-ar wrx e b? -e^ 
K = Eoe —-— Si Sr --■ 

^ph ^ 


( 22 ) 


On closer examination, the long distance behavior 
of gr is mostly determined by the overlap parameters 
Si^r oc. from Eq. (A4), while the short dis¬ 


tance is mostly determined by the hopping parameters. 
However, we note that from a realistic point of view the 
shortest distance in Eig. is given by the size of the 
QDs which is far outside the scope of our approximation, 
oqr > 40 nm for GaAs and ai^r > 18 nm for Si (see Fig.[^ 
in Appendix |B 2\ . Additionally, only a; — yf 0 leads to 
gr yf 0, since, otherwise, the parities of |<i)i) and |<i> 3 ) are 
identical. In Fig. we have plotted the calculated cou¬ 
pling strength gr as a function of the inter-dot distance 
ai for fixed ar/ai = 1.1. Both ai yf and £ 0 are 

required to break orbital and spin symmetries, allowing 
for gr yf 0. We find that a strong qubit-cavity coupling, 
e.g., gr « 2 MHz (red dashed line in Fig. becomes ac¬ 
cessible for an inter-dot distance ai ~ 60 nm, which is 
inside the scope of our approximation. 


IV. IMPLEMENTATION OF LONG DISTANCE 
INTERACTION 


The coupling of the RX qubit to the electromagnetic 
field of a surrounding cavity enables the coherent trans¬ 
fer of information between the qubit system and the cav¬ 
ity. Therefore, if two RX qubits are coupled to the same 
cavity, one can transfer information between them via 
the electromagnetic field. The distance of this transfer 
is limited only by the extension of the cavity. Inserting 
the interaction Hamiltonian Eq. (201 into the two-qubit 
Hamiltonian Eq. 0 , we obtain 


2=1 


Hrw = ^(7^+igi{a+a- a-a^) 


+ Hujpna'a 

(23) 


in the rotating-wave approximation, where we neglect 
the counter-rotating terms due to gr wrx ~ Wph. 
Here a± = ((Ta, ± iCy)/2 are the ladder operators and 
LUi are the resonance frequencies of the i-th RX qubit 
with i G {1,2}. Hence, we can eliminat e the cavity mode 
by a second order SW transformation As a result, 
we obtain 

2 

Htot = Hiq + Hint = 

i=l 

(24) 


\/3 e Eq 

2huJpin 


titrS ^ , 

A 2 — £ 2 “''®! 


dr ) SI Sr 


with the Stark-shifted energy Ceff = huji + 
(21) 5^ (o^a-I-1/2) / [2/i (wi — Wph)] and an effective two- 

qubit coupling parameter g^g = gig 2 {1/ [li(wi — Wph)] + 
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FIG. 6. Doubly logarithmic plot of the gate time tg of the 
photon mediated long-distance interaction as a function of the 
detuning Q — |a;ph — u}i\ The red, blue and green lines 

relate to the correspondingly labeled qubit-cavity coupling 
frequencies gi/{2'Kh) in Fig.|^ For simplicity, we assume that 
the two coupling frequencies fr = gi^ 2 l{ 2 TTh) and resonance 
frequencies a;i ,2 each are identical which leads to tg = 17/8/^. 


1 /[S (a ;2 — Wph)]}, where gi is the coupling strength 
between the i-th RX qubit and the cavity. A review of 
the calculation can be found in Appendix [F} 

The first expression Hiq in Eq. ( [^ yields control over 
rotations around the z-axis of the qubits, while the inter¬ 
action part Hint = ffeff leads to an uni¬ 

versal two-qubit gate. We obtain for the corresponding 
time evolution, 

U(t, to) = exp {-itHint/h) 

1 0 0 0 \ 

0 cos(geffi/^) '^sin (gefft/h) 0 j (25) 
0 isin (gest/h) cos{gest/h) 0 I 
0 0 0 1 / 

which after the time tg = li7r/(2geff) yields the univer¬ 
sal iswAP gate.l^ In summary, we obtain the iswAP-gate 
with the following four steps: 

1. Prepare both RX qubits to be off-resonant from the 
cavity frequency, uji ^ Wph. 

2. Detune both RX qubits to be resonant to the cavity 
frequency, Wi ~ Wph. 

3. Wait for the time tg = h-Kj (2geff)- 

4. Detune both RX qubits to be off-resonant from the 
cavity frequency, uii ^ Wph. 

Together with arbitrary single-qubit gates this allows for 
universal quantum computation over distances on the or¬ 
der of a the extension of the cavity. Fig. shows the 
pulsing time tg, as a function of the detuning = 
|wph —a;i|/27r with i S { 1 , 2 } for specific values of the 
qubit-cavity couplings gi. In order to obtain short gating 
times, we require either a small detuning VLi ^2 or a strong 


qubit-cavity coupling gi. We find gate times tg « 1 ns for 
D « 10 kHz and gi = IMHz which allows for up to 10^ 
gate operations if the qubit decoherence time amounts to 


V. CONCLUSION AND OUTLOOK 

In this paper, we have proposed and analyzed an imple¬ 
mentation of a long-distance coupling of two RX qubits. 
We showed that such a coupling can be achieved by plac¬ 
ing two RX qubits into a high finesse cavity (Q-factor 
exceeds 10 ®) with the condition that the resonance fre¬ 
quency of the cavity matches the energy separation of 
the qubit states. By taking into account the wave func¬ 
tions associated with the RX qubit, we have obtained a 
realistic description of the qubit-cavity interaction which 
yields a microscopic mechanism for coupling the qubit to 
the cavity, namely, by a transition of an electron between 
the outer QDs. In particular, the cavity causes a transi¬ 
tion between states |2) and |3) which are coupled to the 
qubit states by electronic hopping elements. 

For the description of the wave functions and in order 
to estimate the dipole transition matrix element between 
the qubit states we have used the method of orthonor- 
malized Wannier orbitals to construct the electron wave 
functions, by taking only a small, finite overlap between 
the original electron wave functions into consideration. 
As a result, we obtain the qubit-cavity coupling strength 
as a function of the inter-dot distances. Realistic parame¬ 
ter settings yield values on the order of ^ MHz depending 
on the chosen distance. Combining these elements enable 
a coupling between two RX qubits in the same cavity 
with gate times tg on the order of ~ ns depending on 
the detuning between the resonance frequency between 
the cavity and the RX qubits and the qubit-cavity cou¬ 
plings. Realistic parameter choices allow 10® operations 
in the qubit coherence time. 

Since the qubit-cavity coupling mechanism is based on 
the transitions between the asymmetric states | 2 ) and 
13) which are only virtually occupied in the RX regime, 
the asymmetric resonant exchange (ARX) qubilP^ should 
have a strong qubit-cavity coupling, since in this imple¬ 
mentation the asymmetric states are occupied most of 
the time. However, this consideration requires a much 
deeper understanding of the associated wave functions. 
Nevertheless, we encourage the investigation of these as¬ 
pects in future studies, since they may strongly increase 
the qubit-cavity coupling strength. 

This far, we have considered in our analysis only a spin- 
degree of freedom which is appropriate for GaAs. How¬ 
ever, silicon has an additional six-fold degeneracy (two¬ 
fold in typical QDs) of the ground state, the so-called 
valley degree of freedom, due to local inequivalent min¬ 
ima (maxima) in the conduction (valence) band.® This 
leads to a mor e com plex structure of the single-electron 
wave functions . 1 ^ 21111 this paper, we considered a non¬ 
degenerate ground state (strong valley splitting) with the 











same valley ground state in all QDs. In future studies, 
the valley degeneracy can be included, since this addi¬ 
tional degree of freedom co uld p ossibly be helpful, e.g., 
serving as additional qubitsl^^Hl 

In our analysis, we assumed a harmonic confinement 
potential of the QDs. Small anharmonicities can give rise 
to a more complex structure of the wave-functions than 
the one considered here. In future studies, these modified 
wave functions can be included to acquire more accurate 
predictions which will hopefully encourage future experi¬ 
mental implementations of such a setup for long distance 
interaction. 
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Appendix A: Electron wave functions 

Here, we calculate the single electron wave-functions 
of an isolated electron in a QD defined in a 2D electron 
system with a perpendicular external magnetic field B 
applied. We assume that the potential of the quantum 
dot is approximately harmonic. Hence, the correspond¬ 
ing wave functions of the ground state are given by the 
following expression,!^ 


‘p{x,y) 


g-m WQD )/2?i 

ttH 


(Al) 


also obtain a phase shift 


^3{x,y) 


Trh ’ 


^ (SyaTl2l% ^-muiQD [(x-a^)^+y'^]/2h 

Trh 


(A2) 


with the magnetic length Ib = \/h/eB. For the overlaps 
between the left and center QDs, Si = {ipx\^p 2 ), between 
the right and center QDs, Sr = and between 

the left and right QDs, S '13 = we obtain 

Si^r = exp [-a^_^/4as] , (A3) 

S'13 = exp [-(a; + ar)^/4a|] , (A4) 


with the magnetic field dependent QD Bohr radius of the 
electron wave functions 


as 


/ h ^ ma;QD \ 

\4:l^mujQu h ) 


(AS) 


For B = 0, we have og = oq = y^hjmujo. Typical values 
are as « 9.5 nm for GaAs and as ~ 4nm for Si. Without 
loss of generality, we find that the overlaps Si^r and S'13 
are always real due to our choice of the wave functions 
and our assumptions of identical confinement potentials 
wqd in each QD. 

Usually, the overlap S 13 is negligible, e.g., for a sym¬ 
metric setup, ai = Qr, we obtain for the overlaps of the 
neighboring QDs Si = Sr = S with S <C I and for the 
overlap between the left and right QD S 13 = S^ and we 
only have to consider S 13 for very small distances between 
the QDs. 


For B = 0, hujQr) = hujo is the harmonic confine¬ 
ment energy of the QD and typically on the order of 
3 meV for GaAs which corresponds with a QD diameter 
aqD ~ 20 nm, while the higher confinement energy in Si, 
e.g., 6 meV, yields a smaller diameter agp ~ 9nm. Due 
to the magnetic field, B « 310 mT, which is necessary 
to split off leakage state^^, the wave functions are not 
properly described by Eq. (Al) and we have to include 
the influence of the magnetic field. As a result, we obtain 
the Fock-Darwin states, which are the harmonic states 
compressed by a factor b = wqd/wo, where Hujqd is the 
confinement potential of the QD under the influence of 
the magnetic field Here, the modified confinement 
potential is wqd = -^/wq -I- with the Lamor frequency 
lol = eBl2m. Hence, fixing the origin of the coordinate 
system to the center QD, we obtain for the single-electron 
wave-function in the center QD, ip 2 {x,y) = ip{x,y). The 
wave functions for an electron in the left dot are shifted 
by x —>■ a; -I- a; and for the right dot by a; —>■ x — , where 
ai (ar) is the distance from the center dot to the left 
(right) dot. However, due to the gauge transformation of 
the magnetic field A = B{—y,x^ai^r) —> B{—y,x) they 


Appendix B: Orthonormalized Wannier orbitals 

1. Minimizing the localization functional 

To determine the orthonormalized Wannier orbitals in 
the long inter-dot distance approximation we use as an 
ansatz. 


1 ^ 1 ) = (l</5i) +ai |v?2)), 

1 ^ 2 ) = ^(|</52)+a2|v?l)+&2|‘/73)), (BI) 
1 ^ 3 ) = ^ (|</53) +03 Iv?2)) , 

with parameters that are determined by the orthonormal¬ 
ization condition (4>i|$j) = 6ij withz, j S {1,2,3} except 
for one parameter which we define ask = ai Sr/ (a^Si). In 
the next step, we calculate the localization functional^ 
J- = ~ X as a function of k 

which is straight-forward since the Wannier orbitals are 
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FIG. 7. Maximally localized parameter settings for fc as a 
function of the inter-dot distance ai with different dynami¬ 
cally fixed values for ttr- In the symmetric case, ai = ar, 
the optimal setting is fe = 1, while for the asymmetric case, 
ai < Or, the optimal setting is 0.3 < A: < 1 depending on the 
asymmetry and the inter-dot distance. Note that for very long 
inter-dot distances fc —>■ — 1 which corresponds to lihi) —>■ \ipi)- 


superpositions of Gaussian wave functions. The minimiz¬ 
ing conditions, dkJ^ = 0, yields two solutions fc ^ ±1 in 
the parameter regime investigated here, a/ = 45 —120 nm 
(more details are given in Appendix B 2), where we ne¬ 
glect the negative solution, since it corresponds to a van¬ 
ishing overlap between the wave functions. Fig. shows 
the best parameter settings from the minimizing condi¬ 
tion from numerical calculations. In Fig. we plot the 
coefficients of Eq. (B1) as a function of the inter-dot dis¬ 
tance and k = \. 


2. Limitation of our model 

The limitations of the orthonormalized Wannier or¬ 
bitals are illustrated in Fig. |^as a function of the inter¬ 
dot distances. The Wannier orbital wave functions, 
only describe the long inter-dot distance limit. As 
a condition for the validity of our approach, we use 

al^>-2al\og(^--Sl)j. (B2) 

Fig. 1^ shows the resulting limitations for a magnetic field 
of i? = 310 mT for different settings of the confinement 
energy wqd. Considering Swqd = S.lmeV, we obtain 
ai^r > 40 nm in GaAs and ai^r > 18 nm in Si. 


Appendix C: Matrix elements of the position 
operator x 


For an estimation of the qubit-cavity coupling, the ma¬ 
trix elements of the position operator from Eq. (15) in the 


{|0), |1), |2), |3)} basis are needed. Therefore, we express 



ai (nm) 


4>i) 


• 

• 

• 

kl) 

W 2 ) 

ka) 

1^ 

E-2> 

• 

• 

• 

kl) 

Iv^s) 

1^ 

CO 

• 

• 

• 

^ 2 ) 

Iv^s) 


FIG. 8. (a)-(c) Gonstitution of the orthonormalized Wannier 
orbitals. The colored lines illustrate the overlaps Ky’i|4’j)| 
with i,j G {1,2,3} as a function of ai with = l.la; and 
wqd = 3.10 meV in GaAs. 


all three-electron states in terms of our orthonormalized 
Wannier functions. In other words, cj ^ {ci^cr ) creates 
(annihilates) an electron in the z-th Wannier orbital with 
spin a. Hence, we have to calculate the following matrix 
elements, using Eq. (15), 


^) = E 


i,i=l 


(“I \v) 


(Cl) 


with u,v £ (0,1, 2, 3}. For the sake of brevity, we focus 
here on one of the two relevant cases u = 2 and w = 3 for 
the qubit-cavity coupling, 

3 

(2|i|3)= E E Xij (vac| 

i,j = l a=t,l 

Inserting x yields only one non-zero term, since all others 
are zero, since hence, we obtain 

(2|i|3) = -a:i3. (C3) 


The calculations for the other matrix elements can be 
done the same way, in particular, (0| x |1) =0, since two 
electrons have to be transferred (or no electrons have to 
be transferred and a spin has to be flipped). Combing 
all matrix elements, we find the expression in Eq. (16). 
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ai (nm) ai (nm) 


FIG. 9. Limitation of the use of orthonormalized Wannier orbitals in (a) a GaAs TQD and (b) a Si TQD in a magnetic field 
B = 310 mT for different settings of the harmonic potential energies cuqd. Here, the area in red shows the case in which the 
overlap between the QD 1 and QD 3 cannot be neglected and our approximation breaks down. 


Appendix D: Matrix elements of the momentum operator p 


Through the position-momentum relation p = — x] we obtain for the momentum operator in a:-direction 

in the {|0),|1),|2),|3)} basis 




(Dl) 


with the 2x2 blocks 


A = 

B = 

C = 


2V2h 2?i 

^'/^'<T^{'tlX \2 — tlX 21 +t-r{x 23 —X 32 )) 3 im(ti a:i 2 —tj 2:21+tr (2:32— 3:23 ) ) 

2/i 2^/2h 

C \m(tiX21-\/2(A+e)xi2-trXl2-tiX22) -\m(t^X22+tlX2i+\/2(A-e)x22-tj.X22) 

i%/ 3 m{tia:ii —•/2(A+e)a;i2+tra:i3 —*12:22) i%/3m(tra:22 —*12:31-l-\/ 2 (A—£)a :32 — tr2:33 ) 

■ ———————-—5^-—. . 2% " 


i %/2 771*1(2:12 —2:21) 

*1 

im( %/2*,,2:12—4£2 :i 3 — %/2t| 2:23 ) 


im ( V2*,, 2:21-4£2:3 i - G2*1 2:32 ) 

2I 

1^2171*^(2:23—2:32) 

h 


(D2) 


In the case of real-valued matrix elements Xij = Xji € M with i,j € {1,2,3} in Eq. (D2) we obtain a simpler 
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expression for the parameters in Eq. (D1), 


Ar = 



Br = 

Cr = 


im(^tlXii — V2(A+e)xi2 — trXi3—tlX22) i•m(^trX22+tlX3l + ^/2{A — e)x32 — trX33^ 

2 l 2l 

iV3m,(^ttXii — \/2(A+e)xi2+trXl3 — tlX22) iV3m(^tr.X22 — tlX3i + V2{A — e)x32 — tr.X33'j 

2l 2l 


0 

im(^V2t-rXi2—4,5X13 —V2tlX23) 

2l 


im { V 2 t r-2 1 — 4 e a; 3 1 — -\/21 1 £C 3 2 ) 

2?i 

0 


(D3) 


Shifting the momentum operator into the qubit basis of the RX regime acquired through a second-order Schrieffer- 
Wolff transformatiorP^, p = e^pe~^ ^ P — [p, <5'] with the anti-hermitian matrix S. We obtain 


r f A + sB^ -I- Bs^ +B — As -I- sC \ 

P^P-[P,S]- (^5t_st^ + CstC C-stS-^ts j’ 


(D4) 


where s is the 2x2 block in the top-right corner of matrix S. Inserting the expressions for A, B, C, and s we obtain 
the expression p||o),|i) ^ A + sB'^ + Rs'^ given in Eq. (|T7|). 


Appendix E: Calculation of the coupling strength gr 


we have to calculate the matrix elements x. 


In order to obtain the qubit-cavity coupling g, 
i,j G {1,2,3} which we can rewrite due to a transformation back into the non-orthogonal basis of , \(p 2 

^3 


in Eq. (17) with 

|V=3)} 


as Xij = J2ki=i9ki {‘Pk\x\ipi) with k,l € {1,2,3} and qki S K. Due to the separability of the Eock-Darwin wave 


functions in position space, (pk{x,y) = ipk,x{x)Lpk,y{y), where (pk,x{x) G 
are independent. 


the integrations over both space dimensions 


oo / oo \ / oo 

{(fk\x\ipi) = JJ dxdyipl{x,y)xipi{x,y)= | J dx(pk,x{x) x (pi^ 2 ;{x) \ I J dyipl 


,{y)m,y{y) 


(El) 


In the next step we prove that this integral is always real under the assumption that the confinement energy wqd is 
identical in each QD, even for a finite homogenous magnetic field. For the integration over a;, this is true, since the 
wave functions do not acquire a phase in x-direction, hence, it has no imaginary contribution. Proving this statement 
for the second term which is just the overlap in j/-direction, is more complicated due to the additional phase acquired 
by the magnetic field. We obtain, e.g., for A: = 1 and ^ = 2 for the second term, 


dy ^ty{y)v>i,viy) = / dye 


iaiyl2lQ -y^I8as _ 


= dye 


— y^/8a^ 


(E2) 


with y = y — 2iaiag/lg, which is a Gaussian integral shifted along the imaginary axis. Since this integral is real, Xij 
is also real and the only relevant matrix element according to Eq. 0 is 

1 

7Vi N 3 


Xl3 = 


(((^il +ai ((^ 2 ! +bi ((/J3|)i(|(^3) -ba3 \(p2) + h . 


(E3) 
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To calculate ccis, we first consider the more general case of non-vanishing overlap S'13 for our calculation and neglect 
513 later. Inserting the Fock-Darwin functions from Eq. (|A2[), we hnd 


xis = 


mujQD 


dx dy oi exp 


mwQD (a;^ -I- j/^) 


2n 


f mujQu [(a; xf + y^] iaiy 

“n-®^ 


[ mWQD [(flr - x )^ + 2/^] iUry 
+ bl exp <--^:r;75- } I X 


X 03 exp 


2h 

mwQD (a;^ -|- 

m 


211 


+ 63 exp ■; - 


mwQD [(a; -f xf + iaiy 


2h 


21% 


f mwQD - a:)^-I-y^l iary 
- exp <-^ -|- 


2h 


2P 


(E4) 


This integral can be solved straight-forwardly, since all 
integrands are Gaussian. As a result, we obtain 


which yields the expression in Eq. (19) after inserting 
into Eq. ( [I7| ). 


a:i3 = [-az(ai635i -I- aj,Si + 2&3 -|- 513 ) -I- ar{aiSr + 5i3) 
+ biar{cL^Sr + 263 -|- 2 )] / 2 NiN^. 

(E5) 

Setting 5 i 3 = 61 = 63 = 0 and inserting oi, 02, 62 and 03 


from Section III B 1 for the second step, we obtain 


a;i3 = 


Ci\Ci'pSj‘ Qj^dlSl 


2^af -j- 2aiSi + 1-^/a| -I- 203 ^^ -I- 1 
“ arel(nz Ur)5/ 


(E6) 


where dpei can be expanded for k = aiSr/a^Si « 1 (de¬ 
fined and discussed in Appendix |B 1[), 


Orel — 1 + (-5;^ + 5^ -f - ) (fc — 1) (E7) 

Or O/ A , 


Appendix F: Second SW transformation to 
eliminate the cavity mode. 

In this appendix, we present an effective Hamiltonian 
in which the cavity mode is split off by a second SW 
transformation in order to express the cavity mediated 
coupling between two RX qubits.^^ We start with the 
universal Hamiltonian in the rotating frame from Eq. (23) 
and introduce a second order SW transformation 

Htot = -Ho + [iJmt, ^2] /2. (FI) 

where Hq = J2iEoil2 (Tz + hwphO^o and Hint = 
9i i (o’-!-® ~ CT-O^) and the anti-Hermitian operator S 2 
is determined by the SW condition [Ho,52] = —H^t- 
Hence, we express S 2 in terms of operators. 


& = E 


9i 


HuJ-nh — fVjJi 


(cr+a — cr_a^) , 


i=i 

with the previously introduced ladder operators a±. 


(F2) 
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